Stabilizer Formalism for Operator Quantum Error Correction 
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Operator quantum error correction is a recently developed theory that provides a generalized 
framework for active error correction and passive error avoiding schemes. In this paper, we describe 
these codes in the stabilizer formalism of standard quantum error correction theory. This is achieved 
by adding a "gauge" group to the standard stabilizer definition of a code that defines an equivalence 
class between encoded states. Gauge transformations leave the encoded information unchanged; 
their effect is absorbed by virtual gauge qubits that do not carry useful information. We illustrate 
the construction by identifying a gauge symmetry in Shor's 9-qubit code that allows us to remove 
4 of its 8 stabilizer generators, leading to a simpler decoding procedure and a wider class of logical 
operations without affecting its essential properties. This opens the path to possible improvements 
of the error threshold of fault-tolerant quantum computing. 

PACS numbers: 03.67.Pp, 03.67.Hk, 03.67.Lx 



The theory of fault-tolerant quantum computation 
Q 1^ Hi H) IE IS demonstrates the formal possibility of 
efficiently storing and manipulating quantum data for 
arbitrary long times even in the presence of noise, pro- 
vided the noise level is bellow a certain threshold. Fault 
tolerance builds on quantum error correction (QEC) 
P, 0, 0, > which is a mean to actively protect quantum 
information against noise. Encoded quantum states are 
restricted to a code subspace C of the system's Hilbert 
space H = C ® C^. Measurements are performed to 
detect if the noise has taken the system out of C, and 
if required, a transformation is applied to restore it. A 
good code must protect the information against a wide 
range of errors, and admit simple encoding, error correc- 
tion procedures, and fault-tolerant gates. 

Operator quantum error correction (OQEC), recently 
introduced in 0, 1^ , generalizes the standard theory of 
QEC and provides a unified framework for active error 
correction and passive error avoiding techniques such as 
decoherence-free subspaces 0, ^3 ^'^'^ noiseless sub- 
systems jl3, 14, 15J. In this new paradigm, information 
is encoded in a subsystem A of the code space C = A®B, 
and errors need only to be corrected modulo a transfor- 
mation on B. The standard QEC theory corresponds to 
the special case where B is one-dimensional. While this 
generalization does not lead to new families of codes, it 
does allow for new error correction procedures, possibly 
enriching the fault tolerance theory. A prime example 
is Bacon 's OQEC code [H 

that appears to have self- 
correcting properties. 

Most of the QEC codes used for fault tolerance con- 
structions can be described with the stabilizer formalism 
(see Ref. Q and references therein). In particular, the 
first QEC codes proposed by Shor T| and Steane 
are stabilizer codes. Other important examples include 
CSS codes 0,01, topological codes 0, and convolu- 



tional codes [23. A stabilizer formalism has also been 
constructed to describe the passive error avoiding tech- 
niques of decoherence-free subspaces and noiseless sub- 
systems 10. Aditionally, the stabilizer formalism plays 
a central role in other branches of quantum information 
science, e.g. in the so called "one-time" or "cluster state" 
quantum computation model j2]| . Some of the advan- 
tages of the stabilizer formalism are that it provides a 
compact description of QEC codes, admits compact de- 
scription of a restricted class of dynamical systems (the 
so-called Clifford group ) , and allows to build on classi- 
cal coding theory (particularly via the CSS construction). 

In this article, we describe a stabilizer formalism for 
OQEC. We will briefly review the basic theory of OQEC 
and the standard stabilizer formalism. Then, we demon- 
strate a general procedure based on the algebraic ap- 
proach of Ref. [221 to describe the subsystem structure 
A® B using the Pauli group. We also discuss bounds 
that apply to these codes. Finally, we illustrate the sta- 
bilizer formalism by constructing an OQEC code based 
on Shor's 9-qubit code, but which contains a nontrivial 
B subsystem. This code has all the essential features of 
Shor's original code, but admits a simpler error recovery 
procedure and a wider class of encoded operations. 
OQEC theory — Let us first summarize the OQEC the- 
ory. A fixed partition of the system's Hilbert space 
H = A ® B ® is assumed. Information is encoded 
on the A subsystem, i.e. the logical quantum state 
G B{A) is encoded as ® with an arbitrary 

p^ . We say that the physical map £ : B{H) —t B{H) is 
correctable on subsystem A when there exists a phys- 
ical map TZ : B{H) — > B{H) that reverses its action, 
up to a transformation on the B subsystem, i.e. for all 
p"^ and p^ , TZ o S{p^ ® p^) = p"^ ® p'^ for some arbi- 
trary p'^ . In terms of the operator-sum representation 
^{p) — J2a-^apEl, the existence of a recovery map TZ 
requires the following condition to hold (see H, ) 



PElE.P = 1^ (® g^, 



Va, b 



(1) 
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where P is the projector onto the code space — i.e. 



2 



PH = C = A (E) B — and is an arbitrary operator 
in B{B). That this condition is also sufficient for £ to 
be correctable was proven in |23| , along with alternative 
information-theoretic necessary and sufficient conditions. 
As expected, when the B subsystem is one-dimensional, 
Eq. ffl) reduces to the familiar error correction condi- 
tion 50. 

Stabilizer formalism — Let us now focus on the case 
where the system is composed of n qubits, so = . 
The Pauli matrices are defined as 

We denote Xj the matrix X acting on the jth qubit, and 
similarly for Yj and Zj. The Pauli group on n qubits 
Vn is generated under multiplication by the Pauli matri- 
ces acting on each qubits, together with the imaginary 
number i. In terms of independent generators, we have 

Vn = {i, Xi, Zi, . . . , Xn, Zn)- 

The first step in constructing a stabilizer code is to 
choose a set of 2n operators {Xj, Zj}j=i....„ from Vn 
that is isomorphic to the set of single qubit Pauli opera- 
tors {Xj, Zj}j=i^,,,n in the sense that the primed and un- 
primed operators obey the same commutation relations 
among themselves. The operators {X[, ^j}j=i,...n gener- 
ate Vn and behave as single-qubit Pauli operators — we 
can think of them as acting on n virtual qubits. How- 
ever, these virtual qubits have no relation whatsoever 
with the original "bare" qubits (those related to the un- 
primed Pauli operators): the operators X'^ and Z'^ can 
act non-trivially on several bare qubits, i.e. they are 
collective degrees of freedom. It is crucial to note that 
there are several ways of choosing these operators, and 
that these various choices will lead to different codes. It 
should therefore be kept in mind that beside the imposed 
commutation relations, the and Zj are arbitrary. 

The stabilizer group S = {Si,...Ss) with s < n is 
an Abelian subgroup of Vn that does not contain — 1. 
Without loss of generality, we can choose Sj = Zj for 
j — 1, . . . , s. These generators are independent commut- 
ing elements of 'P„, so they can be simultaneously diag- 
onalized. The code space C is the span of the vectors 
fixed by S, i.e. Sj\il)) = |'0) for all j = 1, . . . s, and it 
has dimension 2""". The projector onto the code space 
is denoted P and obviously satisfies SjP = P for all j. 

The normalizer of S, denoted N{S), is the subgroup 
of Vn that commutes with every element of S. (This 
simple definition of the normalizer follows from the fact 
that every pair of element of Vn either commute or anti- 
commute.) Clearly then, elements of N{S) map the code 
subspace to itself. Given the above construction, we see 

thativ(5) = (z,zi,...z;,x^+i,...x;). 

Stabilizer formalism for OQEC — The stabilizer S spec- 
ifies the code subspace C, and we must now define a 
partition of C into subsystems A® B. For this, we fol- 
low the procedure of Ref. and identify subsystems 



via the algebra of operators acting on them.-'^ Central 
to the OQEC theory is a notion of equivalence between 
states: the two states and p'^^p'^ are considered 

to carry the same information even if p^ and p'^ differ. 
To capture this notion, we quotient the code state space 
B{C) by a set of "gauge" transformation Q that defines 
an equivalence relation p ^ p' 3g (z G '■ p = gp'g'^ ■ 
For ~ to define an equivalence relation, Q must have a 
group structure. Clearly, S and i should be in Q as they 
leave states of C invariant under conjugation. For ~ to 
keep states in the code subspace, Q must be a subgroup 
of N{S). Given these properties, Q is a normal subgroup 
of N{S), so C = N{S)/G also has a group structure (the 
quotient group). 

The gauge group Q D {S, (i)} can thus be generated 
by the stabilizer generators, the complex number ?, and 
an arbitrary subset of the Xj and Zj with j > s, i.e. Q = 
{i,Si,..., Sn-r-k, X^^,..., X[^ , Zj^ , . . . , ZjJ where {ik} 
and {jk} are subsets of {s + 1, . . . n). However, for the 
two groups Q and C to induce a subsystem structure on 
C, we must have [Q,C] = (see H^). As a consequence, 
the Xj' and Zj generators of Q must always appear in 
pairs, so without lost of generality, we must have Q = 
{i,Si, . . .S's,X^._^i,Z^+i, . . . ,X^+^,Z^+^) with s + r <n. 
Clearly then, C ~ {X'^+,+^, K+r^i. ■■■,K. K)- With a 
slight abuse of notation, we will henceforth use C to de- 
note the quotient group N{S)/Q and its representation 
on H given above. Since [Q, C] — Q and Q x C ~ N{S), 
it follows from Ref. 22] that these groups induce a sub- 
system structure on the code subspace C = A® B, such 
that the action of any L ^ C and g ^ G restricted to the 
code subspace C is given by 

gP = 11^, ® g^,foT some g^ € B{B), (2) 
LP = L^®]lf,_,, for some L^^ e B{A) (3) 

with A ~ (C^)®*^ and B ~ (C^)®'' as desired. 

To sum up, we have partitioned the n virtual qubits 
defined through the Xj and Zj into 3 sets: s stabi- 
lizer qubits, r gauge qubits, and k logical qubits, with 
s + r + k — n. The Zj operators from the first set are 
denoted Sj with j — I, . . . , s respectively. They are stabi- 
lizer generators and fix the 2'"+'^-dimensional code space 
C. The Z'g^j and X!,_^_j operators from the second set 
are denoted gj and g^ with j ~ 1, . . . , r. They generate 
the group £5 of Pauli operations acting on the r virtual 
qubits of the B subsystem. These qubits do not encode 
useful information: their sole purpose is to absorb trans- 
formations from Q, and as such, they are referred to as 
gauge qubits. Together with the stabilizer and the com- 
plex number i, this set generates the gauge group that 
leaves the encoded information invariant under conjuga- 
tion, Q = Cb X S X (i) . From this definition, it is clear 



While we mostly focus on group structures here, an equivalent 
algebraic description can straightforwardly be obtained by con- 
sidering the associated group algebras. 
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that an Abelian gauge group corresponds to the standard 
stabihzer formahsm, while non-Abelian Q yield OQEC 
codes. Finally the Z'^^^.^^ and X'^^^.^^ operators from 
the third set are denoted Zj and Xj with j = 1, . . . ,k 
respectively. They generate the logical operations C, and 
act only on the k virtual qubits of the A subsystem. 

Although we have given an explicit set of generators for 
C, we stress that only the coset structure of C really mat- 
ters. Operations related by a gauge transformation have 
the same effect on the encoded qubits, e.g any Zj = gZj 
with g ^ G can serve as the logical Z Pauli operator 
acting on the jth encoded qubit. This defines an equiv- 
alence relation Z Z' 4^ ZZ' G G between quantum 
operations. As mentioned in footnote 1, the definitions 
of the gauge group G and logical operations C can be ex- 
tended by considering the associated group algebras: any 
linear combination of elements of G (resp. C) is an oper- 
ator acting solely on the gauge system B (resp. encoded 
qubits A). This extends the notion of equivalence rela- 
tions between states and operations in an obvious way. 

Error correction — We now study the effect of a set of 
errors {Ea} C Vn- Although this may appear restric- 
tive, we note that a recovery procedure TZ that corrects 
{Ea} will also correct any set of errors obtained from 
linear combination of elements of {Ea}. Error detec- 
tion is made by measuring the the stabilizer generators 
Si, . . . , Ss- These give a set of outcomes (mi, . . . , m^) 
taking values ±1, called the error syndrome. The all- 
ones syndrome indicates that the state is in the code 
subspace C, while any other syndrome indicate that an 
error has taken the state out of C. Thus, detectable er- 
rors are those that anti-commute with at least one of the 
stabilizer generator or that have no effect on the encoded 
data. 

To be correctable, the set of errors must satisfy Eq. 
[Note that the conjugation (f) is irrelevant here, so we 
will omit it.] For any pair a, b, the operator EaEi, is an 
element of either Vn — N{S), or N{S) — G, or G- In the 
first case, there exists an S € S for which {EaEi,, S} — 0. 
Plugging in to Eq. ((TJ, we get PEaEbP = PEaEbSP = 
—PSEaEbP — —PEaEbP = 0, so the operator error cor- 
rection condition is fulfilled. In the second case, observe 
that N{S) -G ^ {jC-^I} xG, so Eqs. (EEJ show that 
PEaEbP = L^b ® 9ab fo'" some L^i, ^ U"^, so these errors 
cannot be corrected. For the third case, Eq. (|2J) shows 
that PEaEbP = 11'^ (g) g^f,, so the condition is satisfied. 
Therefore, {Ea} is a correctable set of errors if and only 
if EaEb i N{S) - G for all pairs a, b. 

To construct the recovery procedure, observe that 
equivalent errors Ea ^ Eb have by definition and Eq. (j^J 
PEaEbP = 11^ «) 5fft, and yield the same error syndrome: 
for all 5" G 5 and g G G, [gEa^S] = if and only if 
[Ea,S] — 0. Thus, syndrome measurement can identify 
the coset of {Ea}/G to which the error that occurred be- 
longs. To recover the information encoded in A, we can 
apply any element of that coset to the state. The overall 
effect of this procedure will be a gauge transformation 
since equivalent errors have EaEbP = H"^ <?^b by virtue 
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TABLE I: Stabilizer generators and encoded Pauli's for Shor's 
[[9,1,0,3]] code. 

of Eq. |(5J , leaving the logical qubits A unaffected. 

Bounds — The distance d of a code is given by the mini- 
mal weight^ of operators in N{S) — 5. A code of distance 
d can correct errors on up to {d — l)/2 qubits. A stabi- 
lizer OQEC code therefore has 4 parameters, [[n, fc, r, c?]] 
representing respectively the number of physical qubits, 
the number of encoded logical qubits, the number of 
gauge qubits, and the distance of the code. The Knill- 
Laflamme or quantum Singleton bound n > 2(d — 1) + k 
restricts the possible values of these parameters 0. As 
any bound relating n, fc, and d derived in the context of 
stabilizer QEC, this bound also applies to OQEC. This 
follows straightforwardly from Theorem 3 of Ref. Q . In- 
deed, a [[n, fc, r, rf]] OQEC code can be transformed into 
a [[n, k, 0, d]] QEC code by turning the gauge Z opera- 
tors gj into extra stabilizer generators, i.e. by fixing the 
gauge. 

The theory of OQEC opens the possibility of simpli- 
fying existing codes — turning a [[n, k, 0, (i]]-code into 
a [[n, fc, r, c?]]-code with r > — by identifying "gauge 
symmetries" in their stabilizer. This would lead to more 
efficient error correction procedures with less error syn- 
dromes to measure and wider class of encoded operations. 
A specific example is presented in the next section. The 
key task is thus to find the largest value of r achievable 
given values of n, fc, and d. We have not yet derived a 
general bound for the number of gauge qubits beside the 
trivial observation that at least one stabilizer must be 
measured when d > 0. By exhaustive search however, we 
have ruled out the existence of a "better than perfect" 
quantum code — a 5-qubit code protecting one logical 
qubit against any single qubit error 0, 0] i but which re- 
quires less than 4 stabilizer generators, i.e. that admits 
one gauge qubit. 

Example — Let us illustrate the idea of reducing the 
number of stabilizer generators by identifying gauge sym- 
metries using Shor's [[9,1,0,3]] code [ij. The stabilizer 



^ The weight of an element of Tn is the number of its single-qubit 
Pauli operators in it that differ from the identity. 
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TABLE II: Stabilizer generators, encoded Pauli's, and gener- 
ators of Cb for a [[9, 1, 3, 3]] version of Shor's code. 



generators and encoded Pauli operators for this code are 
given in Tabled 

By inspection, we see that it is possible to combine the 
last 6 stabihzers of this code, thus ehminating 4 of them. 
The remaining 4 stabilizers hence define a 2^ dimensional 
code space, i.e. C contains 5 virtual qubits. However, 
these 5 qubits are not protected against all single qubit 
errors; only one logical qubits is immune to noise, while 
the other 4 extra qubits in C are gauge qubits. Indeed, 
the code defined by Table ITTl is a [[9, 1,4,3]] code, and it 
coincides with the OQEC code defined in JJ]. This new 
code has all the essential features of the original code. In 
particular, it protects one qubit of information against 
any single qubit error, and it has all the features of a 
CSS code (e.g. fault-tolerant transversal c-not). It has 
however lost its ability to protect the logical qubit against 



some 2-qubit errors, but this is not essential to achieve 
fault tolerance by concatenation. Note also that there is 
much more freedom in choosing the encoded operation, 
e.g. the operator II1II2X3II4II5II6X7X8II9 = Sig^X is a 
valid logical X operation. One can easily verify that the 
generators of the gauge group generate Va, so this 
code has 4 gauge qubits as claimed. 

We stress that this code is not the 5-qubits code 0, 0] 
or Steane's 7-qubit code disguised in a 9-qubit code. 
Indeed, the stabilizer of this new code is a subgroup of the 
stabilizer of the original code, and the encoded operations 
are the same as those of the original code. By exhaustive 
search, we have established that the 5-qubit code and 
Steane's 7-qubit code have no gauge symmetry at all. 

Conclusion — Operator quantum error correction theory 
provides a generalized and unified framework for active 
error correction techniques and passive error avoiding 
methods. In this paper, we have developed a stabilizer 
description of such codes. Stabilizer codes have been cen- 
tral to fault-tolerant constructions, as well as other areas 
of quantum information science: it is our hope the gen- 
eralization presented here will enrich these subjects. We 
have demonstrated that bounds which restrict the fami- 
lies of achievable codes derived in the setting of standard 
stabilizer QEC theory apply straightforwardly to OQEC 
codes via gauge fixing. Finally, we have illustrated our 
formalism by identifying "gauge" symmetries in Shor's 
code that lead to substantial simplifications. An impor- 
tant issue which remains open is to bound the number 
of gauge qubits that can be identified given the other 
parameters of the code. 
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[1] P. W. Shor, Phys. Rev. A 52, 2493 (1995). 

[2] A. M. Steane, Phys. Rev. A 54, 4741 (1996). 

[3] D. Gottesman, Ph.D. thesis, California Institute of Tech- 
nology, Pasadena, CA (1997). 

[4] E. Knill and R. Laflamme, Phys. Rev. A 55, 900 (1997). 

[5] E. Knill, R. Laflamme, and W. H. Zurek, Science 279, 
342 (1998). 

[6] J. Preskill, in Introduction to quantum computation, 
edited by T. P. S. H. K. Lo, S. Popescu (World Scien- 
tific, Singapore, 1999), p. 213. 

[7] C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. K. 
Wootters, Phys. Rev. A 54, 3824 (1996). 

[8] D. Kribs, R. Lafiamme, and D. Poulin, Phys. Rev. Lett. 
94, 180501 (2005). 

[9] D. W. Kribs, R. Lafiamme, D. Poulin, and M. Lesosky, 
(2005), quant-ph/0504189. 
[10] L.-M. Duan and G.-C. Guo, Phys. Rev. Lett. 79, 1953 
(1997). 

[11] P. Zanardi and M. Rasetti, Phys. Rev. Lett. 79, 3306 



(1997). 

[12] D. A. Lidar, I. L. Chuang, and K. B. Whaley, Phys. Rev. 

Lett. 81, 2594 (1998). 
[13] E. Knill, R. Lafiamme, and L. Viola, Phys. Rev. Lett. 

84, 2525 (2000). 
[14] P. Zanardi, Phys. Rev. A 63, 12301 (2001). 
[15] J. Kempe, D. Bacon, D. A. Lidar, and K. B. Whaley, 

Phys. Rev. A 63, 42307 (2001). 
[16] D. Bacon, Phys. Rev. A 73, 012340 (2006). 
[17] A. M. Steane, Phys. Rev. Lett. 77, 793 (1996). 
[18] A. R. Calderbank and P. W. Shor, Phys. Rev. A 54, 1098 

(1996). 

[19] A. Y. Kitaev, Ann. Phys. 303, 2 (2003). 
[20] H. OUivier and J.-P. Tillich, Phys. Rev. Lett. 91, 177902 
(2003). 

[21] R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 
5188 (2001). 

[22] P. Zanardi, D. A. Lidar, and S. Lloyd, Phys. Rev. Lett. 
92, 60402 (2004). 



5 



[23] M. A. Nielsen and D. Poulin, (2005), quant-ph/0506069. Phys. Rev. Lett. 77, 198 (1996). 

[24] R. Laflamme, C. Miquel, J.-P. Paz, and W. H. Zurek, 



